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Abstract 

We use an auxiliary field construction to discuss the hard thermal loop effective 
action associated with massless thermal SU(N) QCD interacting with a weak grav- 
itational field. It is demonstrated that the previous attempt to derive this effective 
action has only been partially successful and that it is presently only known to first 
order in the graviton coupling constant. This is still sufficient to enable a calculation 
of a symmetric traceless quark gluon plasma energy momentum tensor. Finally we 
comment on the conserved currents and charges of the derived energy momentum 
tensor. 



1 Introduction 



There has been much discussion on various ways of deriving the energy momentum tensor 
of the high temperature deconfined massless quark gluon plasma (QGP), which is non- 
standard due to the non-local nature of the perturbative resummation required to describe 
the plasma even at the lowest orders of g, the QCD coupling constant 0. Weldon [0 
describes the non-local dynamics of the QGP using a lagrangian with a local dependence 
on auxiliary fields (rather than a non local dependence on physical fields). He then 
proceeds to derive the QGP energy momentum tensor by calculating the canonical energy 
momentum tensor of this lagrangian (which however is neither symmetric nor traceless 
despite the fact that the (classical) field theory of the massless QGP is scale invariant). 
Blaizot and Iancu J| deduce three forms of the QGP energy momentum tensor for a 
pure gluon plasma by integrating the divergence condition d v T^ u = j M (where is a 
current due to an external source). This approach however does not display the link 
between symmetries of the theory and the conserved quantities constructed from the 
energy momentum tensor. 

Brandt, Frenkel and Taylor || discuss the interaction of the QGP with a weak gravi- 
tational field and introduce the two effective actions T HTL [A, g] and Y HTL ^ quark ' ) [A, e], the 
former for gluon graviton interactions and the latter for the interaction of quarks with 
gluons and vierbien fields (where A, g, e are (classical) gluon, metric, and vierbien fields). 
They define the QGP energy momentum tensor to be the sum of the two expressions 

rpHTL = mm ( $F HTL [A,g) \ HTL(quark) = , fT^'lAeK 

which are the definitions of T^ TL and j'^ L \i uark > adopted in this paper. Note that the 
above definition of T^ TL is equivalent to twice the leading temperature contribution of 
the sum over n of thermal Feynman graphs consisting of a graviton current insertion 
T^J -0 -^ and n gluon current insertions. (T^J -0 ^ is the usual gluon contribution to the zero 
temperature QCD energy momentum tensor. It is also the zero temperature coupling to 
the graviton field, tp a ^, which is related to the metric by g a/3 = rj a ^ + K(p al3 where k is the 
graviton coupling constant). A similar interpretation can be given for the quark tensor. 

The properties of the energy momentum tensor such as divergencelessness (unlike 
Blaizot and Iancu we do not have external sources), tracelessness and symmetry (on use 
of the equations of motion if need be) are directly linked to the symmetries of the hard 
thermal loop theory. These properties can be deduced either from the invariances of the 
above effective actions or, equivalently, from the properties of hard thermal loop diagrams 
with graviton insertions. 

In this paper the two effective actions above are written, to first order in the graviton 
coupling constant k <C 1, in terms of local actions with auxiliary fields. This greatly 
simplifies the calculation of T^ L and j^^ L (i uark \ We note that a traceless and symmetric 
energy momentum tensor including a contribution due to the presence of quarks has not 
previously been written to all orders in g in the literature. Furthermore, the auxiliary field 
construction explicity shows that the gluon hard thermal loop effective action (with no 
gravitons), together with the symmetry properties of the gluon graviton effective actions 



given in 0] (and by conditions (A) — > (D) of section (2.2) below), are not sufficient to 
determine the gluon graviton effective action T HTL [A, g], even at O(k). This is contrary 
to assertion made in ||. We give supplementary conditions that enable determination of 
T HTL [A,g] to 0{k). Further, it is shown that both Y HTL [A } g] and T HTL{ - quark \A } g] are 
not known at present to 0(n 2 ). 

We show that the derived gluon contribution to the energy momentum tensor, T B ^ L , 
gives an integrated energy density (i.e. P — j d 3 xT BTL ) which is a positive definite 
functional of the gauge fields, although this property is not manifest. This is a straight- 
forward corollary of the results of Blaizot and Iancu |§ . Finally, it is shown that the gluon 
contribution to the integrated energy density is also given by the Hamiltonian associated 
with the flatspace local auxiliary field action on elimination of the auxiliary fields. 



2 Gluon fields 

2.1 Flat space Gluon auxiliary fields 

The notation used in this paper is as follows: we consider the interaction of a weak external 
graviton field with massless thermal SU(N) QCD in the deconfined phase at temperature 
T and zero chemical potential with Nf quark flavours. The generators of the fundamental 
and adjoint representations of the colour group are denoted respectively by t A and T A . 
These are taken to be hermitian and are normalised by 

tr{t A t B ) = U AB tr{T A T B ) = N5 AB (2.1) 

The structure constants f ABC satisfy [t A ,t B ] = +if ABC t c and (T A ) BC = -if ABC . We 
have quark fields ip, gauge fields in the two above representations A^ = A A t A = A A T A , 
vierbien fields e afl where a denotes a local Lorentz index and graviton fields Q/3 given by 

g afi = r] aP + KLp afi (2.2) 

where k <C 1. 

Initially we deal with the gluon sector of the theory only, although the inclusion of 
quarks will be considered later. We first discuss section (3) of Weldon's paper . Consider 
introducing an action with auxiliary fields Vq(Q,x) = Vq A (Q, x)t A and Wq(Q,x) = 
WQ A (Q,x)t A which transform like the field strength under gauge transformations. Here 
Q is a light-like vector of the form (1, q) and the measure <if2(47r) _1 used below averages 
over all possible q and the label indicates flat space. Weldon defines the action T 
by 

g 2 T 2 X Q {Q- V , Wo, A,rj\ = J ^x2tr [-f m^Wo, + VqQ v F^ v - VgQ ■ DW ,] 

r [V ,W ,A, V ] = g 2 T 2 J^X (Q;V ,W ,A,ri} (2.3) 

where D X W^ = d x W^ + ig[A x , ; m 2 = g 2 T 2 (2N + N f )/18 and £ tr' refers to a 
trace over the colour matrices. X (Q; Vq, Wq, A, rj] is an 'angular integrand' of r with 



(Q; Vo, Wo, A,T)] denoting that X is a function of Q and a functional of the the various 
fields. 

The Euler-Lagrange equations give 



^ = o =► (Q ■ D)Wo fl = Q»F, U ^ = =► (Q ■ D) 2 V , = \m\Q v 1 



We have not yet specified boundary conditions. These are chosen so that the solutions 
we obtain from ( |2.4| ) for the auxiliary fields, when substituted into ( |2.3| ), yield the flat 
space gluon effective action Y BTL [A, rf\. However, we note that T BTL [A, rf\ is not uniquely 
defined in Minkowski space. [As is well known ||, different T BTL [A,r)} correspond to 
different analytic continuations of the uniquely defined Euclidean effective action. This, 
in turn, corresponds to a choice of writing non-local contributions to Tq TL [A, 77] in a 
retarded or advanced form]. Thus, the boundary conditions are not uniquely defined by 
the requirement that we derive a flat space gluon effective action. 

Suppose we use homogeneous retarded conditions (i.e. that Wo, , Vq, — > as t — ► 
—00). Then we must take a retarded solution to ( |2.4|) given by 

Wo, = (Q ■ D); e \(Q»F^) Vo, = (Q ■ D)^ t {Q v F, v ) (2.5) 

where (Q ■ D)~} t is an integral operator whose action on an arbitrary field G A is defined 
by 

((Q ■ D); e \G) A (x) = f° d6U AB (x,x + Q6)G B (x + Q9) (2.6) 

with U{x, x + Q9) AB = (Pexp J° de'[-igQ"A^x + Qe')]) AB 

The resulting expression on substituting (|2.5| ) into (|2.3| ) gives a retarded form of the 
flatspace gluon effective action 

9 2 T 2 x HTL (Q;A, v ] = /rf 4 x2tr[-|mJ((g-D)^g^ A )((g-D) r - 1 t g^ A )] 

rr L [A, V ) = g 2 T*J^X (Q;A,r)} (2.7) 

Alternatively using homogeneous advanced conditions results in a different, advanced, 
expression for Tq TL [A, rj\ above, where the integral operator (Q ■ D)~^ t is replaced in (|2.7| ) 
by (g ■ D)~^ v whose action on an arbitrary field G A is defined by 

POO 

((g • D)~lG) A (x) = - / deU AB (x, x + Q6)G B (x + Q9) (2.8) 



No other boundary condition results in a solution of (|2.4j) , which when substituted into 
( |2.3| ), yields an expression of the form of a possible T BTL [A, 77]. 
We impose a constraint on the gluon fields 

(Q ■ D)£(Q»Fr) = (Q ■ D^KQ-F^) = (Q ■ DY\Q»F, V ) (2.9) 



We denote the set of gluon fields which satisfy ( |2.9| ) by TZq uon . For gluon fields in TZq uon 
the flat space gluon hard thermal loop effective action is uniquely defined and the initial 
choice between advanced and retarded boundary conditions for solutions of ( |2.4p above 
becomes redundant. 

In the following sections we extend the theory to non-zero K and encounter several 
equations for which we have to choose between advanced and retarded boundary condi- 
tions, as for ( |2.4|) above. The different choices lead, for arbitrary fields, to effective actions 
corresponding to different analytic continuations of the uniquely defined Euclidean effec- 
tive action. For simplicity, we impose constraints on the fields such that, as for solutions 
of ( p.4| ) above, any choice between retarded or advanced non-local functionals of the fields 
is actually redundant. 

These constraints imply that whenever in the following sections we encounter a non- 
local expression of the form / d8'B(x + Q6 f ), where B is some function of the fields 
(ip,ip,A,e), the fields are sufficiently constrained to ensure that 

d9'B(x + Qff) = - / d6'B(x + Q6') (2.10) 
-oo Je 

Therefore, we effectively restrict the fields so that in the following sections we can always 
take 

/+oo 
d9' = (2.11) 

This imposes a number of constraints on the collection of fields (ijj, ip, A, e) and we denote 
the set of all such collections by TZ. For example, from above we have that that the gluon 
fields are restricted to some subset of TZ 9 Q luon . The effective actions and energy momentum 
tensors derived below are initially valid only for fields in TZ and are uniquely defined on 
TZ. For fields not in TZ, the effective actions and energy momentum tensors can be derived 
by analytic continuation of the expressions valid in TZ. 

The collections of fields in TZ satisfy a great number of constraints, which are outlined 
below. In the following sections TZ is assumed to be non-empty, unless explicity stated 
otherwise. 

[If TZ is an empty set, the results of this paper are still valid. The calculations, as 
presented below, are invalidated. However, they can easily be repeated but instead of 
assuming the fields are in TZ, we choose some prescription for writing either an advanced 
or retarded form for any non-local expression encountered below (like solutions of (|2.4| ) 
above). Different prescriptions lead to effective actions corresponding to different analytic 
continuations of the uniquely defined Euclidean effective action. The effective action and 
QGP energy momentum tensor derived using this method are exactly the same as a specific 
analytic continuation of the effective action and QGP energy momentum tensor derived 
below. As will be verified below, there is a prescription which yields an effective action 
for fields not in TZ with analyticity properties such that its associated energy momentum 
tensor is retarded. (The physical significance of retarded currents for the QGP plasma is 
discussed by Jackiw and Nair ||)]. 



2.2 Inclusion of weakly coupling graviton field 



Consider a generalisation of the above Minkowskian theory to include interactions of the 
QGP with a weak graviton field, ip a/3 . Following || and || we assume the spacetime to 
be asymptotically Minkowskian. By considering the symmetries of this theory, we aim 
to write down an action of local auxiliary fields which, on elimination of the auxiliary 
fields, gives to 0(k) the generating functional T HTL oi hard thermal loop gluon-graviton 
interactions, i.e. 

T HTL [A, g}= f^Yl • • • ' d 4 x m d 4 Vl , . . . , An (2.12) 

J 47r m,n 

( Xl ), . . . , A^(x m ), <p a ^( yi ), . . . , <p a ^(y n )G HTL (x u /•„.//, y n , Q)Z^a n p n 

= 9 2 T 2 J^X HTL (Q-A,g] (2.13) 

qhtl re f ers ^ Q i\y e h arc l thermal 1-loop truncated diagram with n external graviton fields 
and m external gluon fields. As k < 1 we take n = 1 unless explicity specified otherwise. 
Clearly, the definition of X HTL above only specifies X HTL up to a total angular differential. 
However it is known from hard thermal loop field loop theory that X HTL can be chosen 
to have the following properties 

A. It is non-local but with non-localities only of the form of products of (Q ■ d)^ 1 
operators. Examples of such operators are given by ( |2.3(J| ) 

B. It is homogeneous of degree zero in Q. 

C. It has dimensions of (energy) _2 . 

D. It is invariant under general coordinate, SU(N) gauge , Weyl scaling (and for spinor 
theories) local Lorentz transformations which are restricted to tend to the identity 
at infinity. (A property which all local transformations in this paper are assumed 
to obey). 

When X HTL is referred to in the rest of this paper we implicitly mean the form of X HTL 
which satisfies the above conditions. [An alternative expression differing by a total angular 
differential is the basis of the application of Chern-Simons field theory to hard thermal 
loop phenomena M]. In an attempt to produce X HTL at least to O(k) using an auxiliary 
field lagrangian method, thus enabling a simpler calculation of T^ L , we extend ( |2.3j ) 
to form a general coordinate and Weyl scaling invariant expression. However, we find 
that is not possible to perform this extension uniquely and that this non-uniqueness has 
interesting implications, which will be elaborated below. 

2.2.1 A Weyl scaling and general coordinate extension 

Consider the general coordinate and Weyl scaling invariant extension : 



g 2 T 2 X E (Q;V,W,A,g] = J d 4 x(^g^ e A )2tr 



- lm 2 g W^W v 



+ V^ x F uX -V u [y a V a W» + W a V^y a -y a W a ^] (2.14) 
+E{v u [W x V ll y x - y ■ WdA^ + g„ x W ■ Vy x + W„y ■ OA - y^W a A«}} 
dVL 



T E [V,W,A,g] = 9 2 T 2 J— X E (Q;V,W,A,g] 



where E is an arbitrary constant and Va refers to a space-time and colour adjoint 
derivative i.e. 

V A W^ = d x W„ - T p Xfl W p + tg[A x , W»] 

Note that V acting on a colour trivial object, as will occur below, reduces to the standard 
space time derivative. 

In order to define y fl above, we must consider the null geodesic, denoted by y p (x,9), 
which passes through x, is affmely parameterised by 9 and satisfies the following conditions 

y»(x, 0) = x» y^x, 9) = dy y' °> -> Q p as 9 -> ±oo (2.15) 

d9 

The tangent vector of 0) at 9 = defines the expression ij^ used in ( p. 14 ) i.e. 

(2.16) 



6»=0 



d9 

Note that t/^ generalises Using the null geodesic equation 

y(x, 9) = -!%,{y(x, 0) (2.17) 

and the conditions ( 2.15 ) we can expand ^(x, 9) in powers in of k. Using ( |2.11 ), we have 
to first order in k 

POO 

y x (x,9) = Q x -k d9'-f X (x + Q9') + 0(K 2 ) (2.18) 
Je 

= Q x + k f d9'j x (x + Q9') + 0(k 2 ) 

J —oo 

where 7 A = -Q^Q ■ d)^ + \Q^Q v d x <P^ 

Note that J e ^ 00 d9''~f x (x + Q9') = — d9' , ~f x (x + Q0') gives a constraint on the vierbien 
associated with the metric. This constraint therefore must be satisfied by any vierbien 
belonging to 1Z. Similar restrictions arise at higher orders in k. Without such restrictions 
we can choose y x (x, 9) such that (y x (x, 9) — Q x ) satisfies either homogeneous retarded or 
homogeneous advanced boundary conditions, but not both. 

A is a dimensionless scalar, homogeneous of degree zero in Q, which vanishes when 
g = 7] with the Weyl scaling property that A — > A + 2a. 

[Recall that a Weyl scaling is given by g a ^ — > e 2a g alS , g a p — > e~ 2a g af 3 , A p — ^ A^ , 
if) — > e~%~ip where for this paper we assume a — > as x — > oo] . 



We can give a simple closed form expression for A in terms of y^ix, 9) defined above. 
First recall that a null geodesic, although invariant under a Weyl scaling, is reparame- 
terised so that its affine parameter 9 transforms as dO — > e~ 2 ° ' d9. Thus 

on a Weyl scaling . 

Taking into account the well known (see eg. ||TT|) Weyl scaling properties of the 
Christoffel symbol, it is straightforward to show that 

iy a y a ) - e 2 °(y a y a - 2(y a d a a)) (2.20) 

on Weyl scaling and hence A can be defined by the solution of the equation 

(y-d)A = -(V a y a ) (2.21) 

Again using ( [2.11| ), we have 

A(y(x, 9)) = + d& V u f(y(x, 9')) = - d& V v y»{y{x, 9')) (2.22) 
Je J-oo 

Again the equality of advanced and retarded versions of A requires a constraint on the 
vierbien fields belonging to 1Z. We have 

A(x) = A(y(x, 0)) -> A(y(x, 0)) + 2a{x) = A(x) + 2a(x) (2.23) 

on Weyl scaling as required. It is also straightforward to check this definition is consistent 
with all the other properties required of A. 

Note however A is not uniquely defined at 0(k 2 ) and the consequences of this are dis- 
cussed in Appendix C. Below, we use A as given in ( |2~22l) . However, as we are working at 
O(k) throughout this paper (except briefly in Appendix C), all of the following discussion 
is independent of the second order ambiguity in A. 

Xe{Q] V, W, A, g] above has been constructed to be Weyl scaling invariant if V and 
W are Weyl scaling invariant. [Note that although the Weyl scaling properties imposed 
on W and V appear to be arbitrary , this apparent freedom doesn't actually affect the 
Weyl scaling and general coordinate invariant extension of the flat space action at O(k), 
as demonstrated in Appendix B]. 

To insure Weyl scaling invariance of Xe it is not sufficient to generalise Q-DW^ to y-VW^ 
due to the awkward Weyl scaling properties of the Christoffel symbols present in V and 
hence extra terms are added so that 



Q ■ DW, -> y a V a W, + W a V,y a - y a W a A,^ = y a d a W» + W u d,y» ~ V a W a A 



■A' 



Thus the Christoffel symbol of W u \7^y u cancels the Christoffel symbol of y a V ' a W^ in a 
covariant manner. Using (|2.19| ) and and ( [2.23|) is is straightforward to verify that the 



above expression is simply multiplied by e 2a on Weyl scaling. This ad-hoc fixing of Weyl 
scaling is not unique and an alternative is given by 



Q ■ DW^ = tj^Q ■ DW V -> g^[y a V*W v - W»V^y u - W v y ■ dA + y"W ■ dA] 

The difference between these two possibilities gives the coefficient of E in Q2.14| ) above 
and thus E is an arbitrary parameter, reflecting the ambiguity present in extending Wel- 
don's action. Of course, we don't know yet for what value of E (if any) Q2.14| ) will yield 



X HTL (to O(k) at least) on elimination of the auxiliary fields. Further ambiguities exist 
at 0(k 2 ) and are discussed in Appendix C. 

2.2.2 Elimination of the auxiliary fields 

Let X E (Q; A, g] and r E [^4, g] denote the angular integrand and the action formed by 
the elimination of the auxiliary fields from X E (Q; V, W, A, g] and T E [V, W, A, g] using the 
Euler-Lagrange equations. Thus, 

g 2 T 2 X E {Q-A,g\ = |d 4 a;( v ^^e A )2tr(-|m2^^) (2.24) 

T E [A,g] = g 2 T 2 f^X E (Q;A,g] 

where W is a solution of the Euler Lagrange equations 

y a V a W^ = fF^ + WaiA*)* (2.25) 

(A E )\ = y "A m -v„r + m»y a -y a ^+i9Zy-dA + (a~Li)} 

Clearly any solution to the above equation is covariant under general coordinate 
transformations, transforms like the field strength under gauge transformations and has 
dimensions L . A solution of (|2.25| ) can be obtained as a perturbative expansion in k. 



Writing 

W„ = £ « n W«A* > (y-V) = Q-D+ E« n (2/-V)„, f = Q u + Y, Kn yn ( 2 - 26 ) 

n=0 n ~ 1 n=l 

n=l 

we have 

W , = (Q ■ D)~ l (Q v F IJiV ) (2.27) 

rt-1 
p=0 

Thus Wo p agrees with the flat space auxiliary field in section (2.1). We impose sufficient 
restrictions on the gluon and vierbien fields in 1Z to ensure that the choice of a retarded 
or advanced inverse in fl2.28|) is redundant. 



By expanding (Q-D) 1 in powers of g, it is straightforward to verify that each W n/Ji and 
hence W^ consists of non-localities of the form of products of (Q ■ d)^ 1 . Applying Q • ^ 



to both sides of equation ( |2.28| ) we see (by inducting on n) that each W nil is homogeneous 
of degree zero in Q and hence so is W^. 

W is also Weyl invariant as this is imposed on W when making the extension from 
flat space. Note that any solution of ( |2.28|) must be Weyl invariant as it can be shown, 
by induction on n, that each W nfl is Weyl invariant. 

Hence has the correct properties to ensure X E (Q] A,g] obeys (A) — > (D) as given 
in section (2.2). This fact, together with the observation that in the limit g —> r), X% 
reduces to X , is interesting and contrary to the previous assertion that a functional 
satisfying conditions (A) — > (D) and with the correct flatspace terms would be in fact 
equal to X HTL (Q; A,g] at higher orders in k |6|, This is considered in more detail in 
section (A. 2) of Appendix A. 

2.2.3 The relation between X E (Q; A, g] and X HTL (Q; A, g] 

We show in Appendix A that the following holds: given any expression which 

i. Satisfies conditions (A) — > (D) above. 

ii. Reduces to Xq TL (Q; A, 77] when g = rj. 

hi. Agrees with X HTL (Q; A, g] at 0(g 2 K). 

must, at least to O(k), be equal to X HTL (Q; A,g] for all orders in g. 

Now consider XeLA, g] above. We know it satisfies conditions (i) and (ii) above. We 
will show that for E = , it also agrees with X HTL at 0(g 2 K). First, we calculate S^g 
defined by 

S% = 2g 2 T 2 ^ r , (^S^) 

(so that the angular integral of the above is equal to the energy momentum tensor associ- 
ated with the action Te)- Then we explicitly verify that for E = , S^if is equal to the 
0(g 2 ) contribution to the angular integrand of T^ g TL which is consistent with conditions 
(A) — > (D) (as calculated using thermal field theory in ||). This shows equality between 
Xe and X HTL at 0(g 2 n) which explicitly confirms condition (iii). Hence we deduce that 
X E=0 (Q; A,g] = X HTL (Q; A,g] at least to O(k) and hence that J dVL(Any l Sfj is equal 
to the gluon contribution of the retarded QGP energy momentum tensor, T^ g TL , to all 
orders in g. 



2.3 Calculation of S^g and the QGP energy momentum tensor 

The calculation of S®a is simplified by noting that for when V, W are considered as 
functional of the physical fields ie. V = V(Q; A,g], W = W(Q; A, g] we have 



SX E (Q;V,W,A,g] 
5^ 



Ayw 



5X E {Q;V,W,A,g] 
8X E (Q;V,W,A,9] 

SWn 



SX E (Q;V,W,A,9] 
SW^Q\A,g] 



SV^Q;A,g] 

A,g,W 



A,g,V 



6g° 



but with g = n, V = Vq, and W = W (see 0) we have sx e(Q;V,w,a,9] and gx B (Q|Vjy,A, g ] 
zero and hence 

, XYnfO-V = VJA a\ W = WJA n] A n] 

(2.29) 



A,V,W const 



Note that calculating S®a only requires the flat space equations for the auxiliary fields and 
thus that Xe{Q] A,g] at O(k) depends only on the flat space auxiliary fields. Performing 
the calculation it is useful to note the following 



d 4 xG(x)((Q ■ d); e )F)(x) = - J d*x((Q ■ d)-lG)(x)F(x) (2.30) 
d9F(x + Q6) ((Q • d)~} v F)(x) = - / d6F(x + Q6) 

-oo JO 

The labels 'ret' and L adv' can be dropped if which / d8'F(x + Q6') and / d8'G(x + Q8') are 
non-local functionals of the fields for which ( 2.11|) holds. We impose sufficient restrictions 



for fields in 1Z to ensure that all non-local functionals encountered below satisfy (|2.11 ). 
Hence from the definitions of y and A, we have for any F\ and F 

Jdxy x F x = JdxrmQiaVpy-QaQpiQ-d)'^^ (2.31) 

J dxAF = J dxr P \[Q a Q^(Q ■ d)- 2 + Vaf) - Q {a dp } (Q ■ dy^F 
where A{ a Bp} = A a Bp + A B a 

On performing the calculation we find 

S$ = 4tr(|[Q a Q /3 D(Q • d)- 2 - Q {a d n (Q ■ ^-^[-fmJW^WoJ - \m 2 g W^ a W^ 

+i[Q{aVpy - QaQp(Q ■ d)- l d x }[V£F^ - V Qn D x W^ + d^WuYg)] (2.32) 
+\e{[Q ■ d{W m V 0a} ) - [Q {qV x } - Q a Q p (Q ■ d)- l d x \[d^W QX vVg + V QX W$)\ 

-[QaQMQ ■ d)- 2 + - Q {a d n (Q ■ dy l ][Q ■ d(y 0lt wg)]}) 



Note Vqu and Wo u are given by fl2.4|) and we have used Q^Vq^ = Q^Wq u = in deriving the 
above. The the coefficient of E in ( 2.32Q is discussed further in section A. 2 of Appendix 



Setting E=0 we find gives equation (6) of in momentum space at 0(g 2 ) thus 
showing, given the result of Appendix A, that Xe=o(Q; A, g] = X HTL (Q; A, g] at O(k) for 
all orders in g. Thus is the angular integrand of the gluon contribution to the QGP 



energy momentum tensor T^ L and hence 

rpHTL _ [ dVt qE=0 , < 

af3 ~ J 4vr a/3 1 ' 

The calculations deriving this are valid for gluon fields in 1Z. We can analytically con- 
tinue from 7Z to consider this expression for gluon fields not in 7Z. A choice of analytic 



continuation where all fields are retarded is of particular physical significance ||. Despite 
using a local auxiliary field angular integrand, T c ^ Ti is not local even with respect to the 
auxiliary fields. This is because ?/ M and A have a non-local dependence on the graviton 
fields. 

An interesting but difficult question is whether or not Xe=o{Q', A, g] is correct at 0(k 2 ) 
and some of the difficulties that arise are discussed in Appendix C. 

[We again consider the problem of choosing a prescription which yields a retarded 
energy momentum tensor. As discussed in section (2.1), such an approach is required if 
we drop the assumption that TZ is non-empty. The prescription is to choose the advanced 
expressions for y and A in ( [2.18 ) and ( 2.22 ), while choosing retarded expressions for all 



other possibilities. With this choice of prescription we see that, on use of ( |2.30|) , ( [2.31| ) 



is still valid on replacing (Q ■ d) 1 by (Q ■ d) r ^ t . This would therefore lead to retarded 
energy momentum tensor, automatically valid for fields not in TZ}. 

3 Single auxiliary field model 
3.1 Flat space auxiliary fields 

The same analysis is applied to a model with only one auxiliary field m{J (Q, x) = itiq A (Q, x)t A 
which again transforms like the field strength on gauge transformations. 

Consider the following angular integrand X (Q; m , A, rj] and action T [m ,A, rj\, as 
flat space precursors: 

g 2 T 2 X (Q;m ,A, V ] = J d 4 x(4b)tr[m*Q v F Xu + \{Q ■ Dm 0u )(Q ■ Dm")} (3.34) 
T [m ,A,7]] = g 2 T 2 J ^X (<2;m , A } rf\ 



On taking b = \m 2 and using (|2.11|) together with the Euler-Lagrange equation for m, the 



above reduce to the non-local flat space hard thermal loop angular integrand and action 
i.e. 

sx 

5m 

g 2 T 2 X HTL {Q . Av] = _ Jd * x 2tT[lm 2 [(Q.D)-^F,MQ-D)- l Q x F\}} 



m 0tl = (Q-D)- 2 Q»F^ (3.35) 
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Note that above we have used 



T [A,rj}=g 2 T 2 S^X (Q;A V } 



J d 4 x(Q u F fl v )[(Q ■ DY 2 (Q V F^ U )} = - J d 4 x[(Q ■ D)-\QT , F ltt/ )][(Q ■ D^iQ" 'F\)|3.36) 
which is valid for fields in TZ. 

3.2 Inclusion of Weakly coupling graviton field 

Again, on generalising to curved space and imposing Weyl invariance , we hope to obtain 
X HTL [A, g] on the elimination of the auxiliary fields. As in section 2 we find there are 



ambiguities. Consider the generalisation where E' is an arbitrary constant: 



g 2 T 2 X E , = J d 4 x(^ge K )2ti\2b(T] u p mPy x F vX ) (3.37) 

+b g^{B p B v + E'B^C" + ±E' 2 C»C U ) 

where B^ = (y a V a m p - m x V \ij p + y p m x d x A) (3.38) 
C = (m p V u y p - m ■ yA, p +m u y ■ VA + m ■ Vy u - y u m ■ VA) (3.39) 

with m p — > m p on Weyl scaling. Again the choice of the Weyl scaling properties doesn't 
affect Xe'(Q', tu, A, g] at O(k) (see Appendix B) and the choice above is for simplicity. 

The coefficients of E' & E' 2 correspond to to the fact that the extension to a general 
coordinate and Weyl scaling invariant action is not unique and they are due to taking 
different combinations of the following two possibilities for the generalisation of Q ■ Dm} 1 
(both of which are simply multiplied by e 2cr on Weyl scaling). 

Q-Drn^ -> ij a V a m p - m a V a y p + y^m a d a A 
Q • Dm? = tfQ ■ Dm v -> g^[y a V a m u + m p V u y p - m a y a A,„ +m u y a d a K\ 

The Euler-Lagrange equations for m^ are 

y ■ V(B a + E'C a ) = y x F aX + (B u + E'C u )(/\ E 'Y p (3.40) 

where (A E ') p is given by ( |2.25|) (for E — ► E'). Thus by ( |2.25[ ) we see that W p solves 
(|3.40|) and hence we have 

y a V a m p = g^W v -{A E 'Y v m v (3.41) 

This is solved in exactly the same way as ( |2.25| ) (after raising indices and replacing y v F ,J ' v 
by g pl/ W v ). Clearly any solution m M to the above equation is covariant under general 
coordinate transformations, transforms like the field strength under gauge transformations 
and has dimensions L°. By expanding in powers of g and k it is straightforward to 
show that it non- localities are only of the form of products of (Q ■ d)^ 1 . Writing m M = 
£ n=0 K n m^,it is also straightforward to prove [by induction on n exactly as in section 
(2.2.2)] that m^ is homogeneous in Q of degree minus one and that its imposed Weyl 
invariance is consistent with the Euler Lagrange equations. 

Defining X^i(Q; A, g] by (|3.37|) at the solution of t \3.40j ) we see that m M has the correct 
properties to ensure that X™,(Q;A,g\ satisfies condition (i). Section (3.1) shows that it 
satisfies condition (ii). Thus we see that X^,(Q; A, g] is a counterexample to the assertion 
that a functional satisfying conditions (A) — > (D) and with the correct flatspace terms 
would be in fact equal to X HTL (Q; A,g] at higher orders in k || [7|. [In fact, as we will 
show below Xf,(Q; A, g] = X E (Q; A, g] at O(k) for E = E'}. 

In addition condition (iii) is satisfied for E' = 0. We show this by simply calculating 

*»S* = 2g 2 T 2 Ife^ ggpdld (3 .42) 



and verifying that m S^p =0 = S^g° (which is equal to the angular integrand of T^p L which 
satisfies (A) — > (D) by the results of section 2). 

Using the techniques of section (2.3) we find m S^ /3 is a function of the gluon fields 
given by 

m S% = Abti{\[Q a Q p U(Q.d)- 2 -Q {a d p] {Q-d)- l }[m ^ + \Q-Dm^Q-Dm^ 

+ \[Q{cWp} - QaQp(Q ■ d)- l d x ][m^F\ + (D x m 0fl )Q ■ Dm% (3.43) 

+d> x (m 0lJl Q ■ Dm ox )} - \{Q • Dm 0a Q ■ Dm 0/3 ) 
+ \E'[- [Q {a 7i$ } - Q a Qp(Q ■ d)- x dx\[dp{Q ■ m ox m% + m 0X Q • m£)] 
- [QaQMQ ■ 9)- 2 + Vap - Q{«d p} (Q ■ d)- l }[Q ■ d(m Qfl Q ■ <)] 
+ Q ■ d(Q ■ m 0{f3 m 0a} ) } 

where /** = Q P F^ p and m ^ is given by ( |3.40| ). Note that we have used Q^m^^ = 
in deriving the above and that the coefficient of E' in m S^p above is the same as the 
coefficient of E in ( |2.32j ) which is consistent with uniqueness results stated in section (2) 
(and derived in Appendix A). 

Setting E' = in ( |3.43| ), applying the identity 

(D x m 0fl )Q ■ Dm 0lM = d x {(Q ■ d)-\Q ■ d)[(Q ■ Dm%)m Qtl ]} - m%D x (Q ■ Dm Q J 

in the second line above and noting that Wo M = (Q ■ D)m 0fJi and Vo M = 2b Q ■ Dm 0fJi it is 
straightforward to see that m S^ /3 =0 is equal to S^p derived earlier. 

Hence X^, =Q (Q; A, g] satisfies conditions (i) — > (iii) . Indeed, we see from the equiv- 
alence of S^p and m S^ f3 (for E = E') that the actions presented in section (2) and this 
section are equivalent at O(n) and hence that the results of section (2) could be derived, 
ab initio, with a model using only one auxiliary field. 

Again, considering agreement of these actions and the gluon graviton hard thermal 
loop effective action at 0(k 2 ) is a more difficult problem. 



4 Quark Sector 

4.1 Flat space auxiliary fields 

Again we perform the analysis of section 2 only now for the quark sector of the theory, 
with the aim of deriving an expression for 

a 5T HTL ^ uark ^[ilj,ilj,A,g]' 



linie^ 



e 



5e af3 



{AAA) 



where e is the vierbien field. 

We start from a flat space angular integrand X quark ^ and action T quark ^ with a 
dependence on the auxiliary fields Xo(Q, x ) an d *&o(Q,x) , the latter of which acts as a 
lagrange multiplier (cf. V(Q,x) in section 2). 

g 2 T 2 X quark{0) = J d A xa(ipxo + ^o(Qip + iQ • Dxo) + h.conj.) 

r — ^2^2 r dQ y- 

1 quark(0) ~ 9 1 J -^^quark(0) 



where a = g 2 T 2 (N 2 — l)/(16iV) and D x is the covariant derivative acting on fields in the 
fundamental representation ie. D x ip = d x + igA x ip. The Euler-Lagrange equations give 

(Q • D) X o = i<%i/> and (Q ■ D)^ = iip (4.45) 

Again, we impose sufficient restrictions on the fields belonging to TZ to ensure that ( |2.f f| ) 
holds for the non-local functionals of the fields which we encounter in the quark sector. 
Thus we have 

Xo = i{Q • D)- 1 ^ = i(Q ■ (4.46) 

with {Q-D)-^ = (Q-D)-^ = (g-D)-V (4.47) 

where (Q ■ D)~^ t and (Q ■ D)alv are integral operators similar to (|2.6|) and given by 



((Q ■ D);^) A (x) = f° d6U(x,x + Q9)^(x + Q9) (4.48) 

roo 

((Q ■ D)-} v i,)\x) = - d6U(x,x + Q6)tlj(x + Q6) 

Jo 

with U(x,x + Q6) = (Pexp f d6'[-igQ^A^x + Q6')}) 

J 

Thus we have on elimination of the auxiliary fields and using [Q, Q ■ D] = 0, 

X i0)guark = ia$Q{Q ■ L>)~V + h. conj. (4.49) 

This is the quark flatspace hard thermal loop angular integrand denoted X^^ k ^ and the 
above calculation is valid for fields in TZ. We can consider (|4.49| ) for fields not in TZ by 
analytic continuation. 

4.2 Inclusion of weakly coupling graviton field to quark sector 

Again, we consider generalising to curved space and imposing Weyl invariance noting 
that the Weyl transformation properties of the auxiliary fields once more do not affect 
the action at O(k) (see Appendix B). Thus \1/ and x can be considered, without loss 
of generality, as Weyl invariants. Unlike sections 2 and 3 we find no ambiguities on 
generalising the angular integrand and action. This is consistent with the fact that quark 
sector analogues of the gluon conditions (i) and (ii) (given by (i)' and (ii)' in section (A. 3) 
of Appendix A) are sufficient, without an analogue of the gluon condition (iii) above, to 
determine x HTL( - quark \ (See section A. 3 of Appendix A for further details). This makes 
the quark analysis relatively straightforward. Y HTL ( quark ^ and x HTL ( quark ^ are defined 
analogously to their gluon counterparts ie. 

vHTL (quark) = f ^^4^^ ^ _ ^ d A x m d A y x , . . . , d A y n A w (xi ) , . . . , A^ m (x m ) 

J m,n 

^(yi), • • • , ^{ Vn mz x )G HTL ^ rk \x u r„. !h y n] Q)^ anP J{z 2 ) 

= g 2 T 2 J^X HTL ^ k \Q^^,A,g] (4.50) 



where Q HTL (i uark ) refers to the hard thermal 1-loop truncated diagram with a quark anti- 
quark pair of external fields, n external graviton fields and m external gluon fields. As 
before we take n = 1 and assume that x HTL( - quark ^ satisfies conditions (A) — > (D) given 
above. Consider the Weyl scaling and general coordinate invariant generalisation of the 
auxiliary field angular integrand in section (4.1) : 



g 2 T 2 X quark = Jd 4 x^a(ije r fx + *(e^(e"j + ie A y ■ V X ) + h. conj. (4.51) 

r quark J ,i_ quark 

where 



(y-V)x = y v (d v + igA u + J ubc a bc ) X (4.52) 
with u' ubc = +\e p b e cp . v + \{e bu e p c - e cu e p b )A, p (4.53) 

and e~$i = {y v e av ~j a ) (4.54) 

Note that e A (y ■ V)x and e~%$ are both Weyl invariants. 

We now eliminate the auxiliary fields x, \I/ from (|4.51 ) and denote the resulting ex- 
pression by X* uark . We have from the Euler-Lagrange equation for \I/ 

^fe = ^(y■V)x = ^[e^e-^)^ (4.55) 

Again, we restrict the fields so that we have equality between the retarded and advanced 
solutions of fl4.55p . Thus, 

/0 p oo 

V(x, x + Q6) x (x + Q6)= V(x, x + Q6) X (x + Q6) 
-oo JO 

where V(x, x + Q6) = (P exp - f d6'[y l '(igA v (y(x, 6')) + u' ubc (y(x, 6'))a bc )})(4.56) 

J — Q 

= (Pexp- fjigyMvW)))] 

(pexp- [yyuUM^^yi 

where y(x,9') is given by (|2.17|) subject to ( |2.15|) . 

Noting that \l/ and \I/ are Lagrange multipliers, we can now eliminate the auxiliary fields. 
Using [e~2^ , (j) ■ V) -1 e = for any \ (See Appendix D), we have the following 
expression for X* uark 

t, = ./A»*ft.Vr'^4b. con, ) (4.57) 

It is straightforward to verify that X* uark satisfies conditions (i)' and (ii)' of section (A. 3) 
of Appendix A, where it is shown that this implies that X* uark = x HTL( - quark ^ to O(k). 



Therefore, we can use the same methods as before (section (2.3)) to calculate the quark 
energy momentum tensor : 



t: 



HTL(quark) 



a/3 



J ^Sl u f k (4.58) 



xv* A Y HTL ( ouark ) 

wither - 9 2 T> lim^fe-^*) =g*I* Bin^fe ™ ^ 

S q a J rk = alQaQ^iQ-dy'-Q^iQ-d)- 1 }^^} 

+ a[Q {a r) X p } - QaQpiQ ■ dyW^oj^ + i* Dx<8* ] 
+ ia[d^Q a ap M x + ^Qp°a M X + ^Q^apx)] + h.conj. 

It is straightforward verify that, up to O(k), X* uark agrees with the quark gluon 
graviton effective action given by equation (34) of ||. They do not agree at 0(n 2 ) despite 
the fact that both satisfy conditions (i)' and (ii)' due to the fact that our definitions of A 
are different at 0(k 2 ) (as discussed further in Appendix C). 

Again, these calculations are valid in TZ and may be analytically continued to fields 

off n. 



5 Conserved Currents 



First of all we note that for any G\ 

J d 3 xG = J d 3 x{Q ■ d)~ 1 d x G > 



(5.59) 



Using ( |5.59| ) it is straightforward to calculate the integrated energy and momentum den- 
sities. Consider the gluon sector. We have from fl2.32|) (with E = 0) 



d 3 xT t 



HTL 
Oa 



dQ ,o , 
— efa;4tr 

47T 



h*°V»» - ^ou){\m 2 g W^) (5.60) 



where Wq and Vq are given by ( |2.4|) . Notice that P a is local before the elimination of 
the auxiliary fields. It is straightforward to verify that this expression agrees with P a as 
would be calculated in other papers (compare eg. (3.13) of 0) which Blaizot and Iancu 
[0], using properties of the angular integral / dfi(47r) _1 , have shown to be equal to 



Po = / ^x \m 2 g [{Q ■ Dy^FoMQ • D^Q^} 



(5.61) 



Thus Po is i n fact a positive definite functional of the gluon fields even though this 
positivity is not manifest. 

Similarly, we can write down pv uark 



pquark 



d 3 xT HTL( - 9uark ^ 



(5.62) 



dfl 1 ; 

— I d xa 

47T 



Vaoi^oQ^) + ^o7«^ + iV D a QV + h.conj. 



where \l/o is given by (|4.45|) . We note that this expression is a local function of the auxiliary 
fields and does not appear to be positive definite which perhaps is to be expected eg. by 
comparison with the zero temperature quark contribution to energy momentum tensor. 
T^ L and j^ TL ('?" arfc ) are symmetric and traceless (on use of the equations of motion 
for T HTL<yquark ' > ) (see |10|). Thus, the construction of angular momentum and scaling 



currents/ charges is standard. They are defined by (for the gluon sector) 

yf , — rn ( T HTL _ rpHTL\ ft _ 1/rpHTL 

Mn\v — X^[l Xu x ^ 1 Xii I A* 1 fiv 

J„ v = [ d 3 xM^ v C= I d 3 xC (5.63) 



with similar definitions for the quark sector expressions. It should be noted that these 
expressions unlike the P M are non local even before the elimination of the auxiliary fields. 

As a final exercise, we compare Pq = J d 3 xT^ TL with the the Hamiltonian derived 
from the local auxiliary field lagrangian density C where g 2 T 2 Xo(Q; mo, A, rj\ = J d A xC. 
We use the m{J field model for simplicity as it has less constraints then the Vq, Wq model 
where V$ acts lagrange multiplier. 

In the following we write we write out colour indices (A, B, . . .) instead of using a 
matrix notation. For clarity, we also drop the flat space label o although we are working 
with flat space fields in this analysis. Thus all indices in the remainder of this section are 
either colour or coordinate ones. 

£ = (2b)[m XA Q^F A + \{Q ■ Dm A )(Q ■ Dm vA )\ - \F A F^ A (5.64) 

The conjugate momemta for m and A are given respectively by 

W A = 2b(Q ■ D)m^ A W A uon = E^ A + 2b{Q»m 0A - m» A ) (5.65) 

where E^ A = F^ A . Hence 

d m^ A = ^W A + Q-"Dm^ A + if ABC A B m^ c (5.66) 

d Af = DjA A — E A n;f uon = 

U = U A d m^ A + U^doAi - C 

= |(E 2 + B 2 ) + U A (W A + Q ■ Dm" 4 + if ABC A B m» c ) 

- 2bm' A Q k F A + A A (D j U A ) (5.67) 

We want to compare / f- J d 3 xH with f d 3 x m S^=°. 



4-7T 

[Note that in deriving m S E f 3 =0 in (|3.43|) we have used the equations of motion for 



m 



m qE'=0\ 



to drop a term with a (Q^m^) dependence, which we denote by S( m S^g j 

8{ m S%= Q ) = -llQ^^Q-dy'+v^-Q^yiQ-d)- 1 } 

Mm^iQ ■ Dm p ))] 

However it is straightforward to check that J d 3 xS( m SQQ =0 ) = and hence / d 3 x m S^ =0 is 
given by integrating ( |3.43| ) with E' — 0]. 



Subtracting / *p / d 3 xH from J d 3 x m 5' ( ^ ) ' we obtain not zero but 

J^J d 3 x2b m 0A [(Q ■ D) 2 m 0A - Q»F 0A ] - A^DjUf + if BAC U*m» c ) (5.68) 

Thus the m field theory produces a non standard result due to the unusual form of its 
'kinetic energy' term |(Q ■ Dm A )(Q ■ Dm uA ) which has cross terms of time and and spatial 
derivatives. However, the physical theory is obtained by elimination of the auxiliary fields 
in which case the term proportional to m° drops out (see |3.35| ) leaving only the Gauss 



constraint term, ie. a gauge transformation generator which is constrained to zero when 
considering the physical, gauge invariant, theory. Hence, we see that for the physical 
theory J § f d 3 xH = f d 3 xT* TL . 



6 Conclusions 

In the above we have shown how to derive the gluon and quark contributions, T^ L 
and r p^ L ^ uark ^ t the retarded QGP energy momentum tensor using an auxiliary field 
method. These tensors have some desirable properties : 

a. They are automatically gauge invariant, without having to add divergences 'by 
hand'. 

b. They are symmetric. T^ff L = T^J L from its definition and j^ TL ('? MarA: ) _ rpHTL(quark) 
as a consequence of the local lorentz symmetry and the equations of motion (See 
eg- E 



c rpHTL ^ a _ rpHTL{quark) ^ a _ g us j n g w e yi scaling symmetry together with the 
equations of motion. 

d. d a (T^ L + j<HTL(quark)^ _ ^ &g required for a energy momentum tensor of a closed 
system. This is a consequence of general coordinate symmetry and the equations 
of motion. (Note that we have not introduced an external source in this paper and 
hence do not a external source current on the right hand side of the above, unlike 
Blaizot and Iancu). 

e. Despite appearances, T^ L is a positive definite functional of the gluon fields, al- 
though to show this requires properties of the angular integral ie. / dQ^ii)^ 1 . 

Properties (a) — > (c) allow straightforward construction of the conserved quantities of the 
QGP ie. momentum, angular momentum and the scaling current as given above. 

However, the expressions for T^ L and r p^ L ^ uark ) are cumbersome, being non-local 
functionals even before elimination of the auxiliary fields. The angular momenta and 
scaling currents/charges are non local even before elimination of the auxiliary fields mak- 
ing their manipulation cumbersome too. Finally we conclude by noting that the insights 
gained using the above method have shown that T HTL [A, g] and T HTL ^ quark ^[tjj,ip, A, g] 
are more complicated than previously thought and are certainly not known at present to 
0(k 2 ) . 



A Appendix A 



In this Appendix we discuss the result that any expression, denoted X 9uess (Q; A, g] which 

i. Satisfies conditions (A) — ► (D) 

ii. Reduces to Xq TL (Q; A, rj\ when g = rj 

iii. Agrees with X HTL (Q; A,g] at 0{g 2 n) 

must in fact equal to X HTL (Q; A, g] at O(k) for all orders in g. In appendix C we briefly 
discuss the some of the further complications that arise at 0(k 2 ). 

Firstly, we note that this result is much weaker than previous speculation that (i) & (ii) 
would be sufficient to imply equality between X 9uess (Q; A, g] and X HTL (Q; A, g] [|, 0]. 
As seen at the end of section (2.2.3) Xe is a clear counterexample. 

As condition (ii) gives equality at 0(k°) we are primarily interested in the theory at 
O(k). It is sufficient to consider 

SiJ ss = 2^T 2 lim^ ( 5xau l$^ ) (A.Q9) 

Suppose two different X 9uess (Q; A, g] exist satisfying (i) — ^ (iii) above, leading to two dif- 
ferent S 9 ^f ss whose difference we denote by AS a p- 

We show below that (i)— ►(iii) force AS a p to zero at all orders in g which proves the result 
stated at the start of this appendix. 



A.l Conditions (i)^(iii) =^> AS a/ 3=0 

Firstly note that (iii) implies that AS a/ 3 is zero at 0(g 2 ). The idea of the proof at higher 
orders in g is simple and in the same spirit as 0. Write down the most general 0(g s ) [here 
s > 3] contribution to AS a p consistent with conditions (A)— ►(C) of (i) above. Next, in a 
step that requires a lot of tedious but not difficult algebra, show that the homogeneous 
Ward identities that correspond to condition (D) of (i) above can only be satisfied if AS a p 
at this order is identically zero. 

The most general 0(g s ) contribution to AS a/ 3, denoted AS? , [where a — ► /x s+ i 
and (3 — ► fi s +2 ] can be written as the following integral 

A^ s+lMs+2 = Jd%...d%d i Pl ...dW^ Pi<yi - x) (A.70) 

K,..^s + ^ +2 Afll (yi)--- A " s (ys) 

where Pi (i — 1 . . . s) denote the i th external gluon momentum. W s is a function of pi . . . p s 
given by the expression below. First, some definitions (used only in this appendix). 

Qcf^r- ~ Qpi • • • Qui-iQui+i • • • Qnj-iQuj+i ■ ■ ■ Qfj.s+2 

Hifij... 

A ijkl = ^ijklVuiUjUkUiQ^—^—-^ Brf = Bijll^-^.Q^y^ 

P'iH'j 



where A?fcz, -By, (Cy'V* are (thusfar arbitrary) coefficients constructed from the gluon 
momenta and their contractions with Q M and each other. Note that the index of 
{Cij)n k corresponds to an index on some combination of the external momenta (and thus 
it does not belong to a or 77). D^. and E^. below are defined in the same way. Thus 
W s can be written in form compatible with (A) — > (C) as 



• • -Hs+2 — {(^1234 + ^1324 + ^1423) + (^1235 + ^1325 + ^1523) + •• • (A-71) 

+ (As-l,s,s+l,s+2 + As-l,s+l,s,s+2 + As,s+2,s-l,s+13) 
+ B* 2 + . . . + -B* +ljS+2 + (Ci 23 + C 231 + Cg 21 ) 

+ . . . + (C s _i s s +i + C SjS+1jS _ 1 + Cg+i s s _i) + -D AllAt2 Q / _^_ N 

/il/i 2 

+ • • • + D^ s+1 ^ s+2 Q, » v + E^Q^^ + . . . + £? /is+2 <5/-^-s| 

The above expression is over-general in that it does not take into account all the 
symmetries of W s , such as symmetry on interchange of graviton indices (/t s +i «-> /is+2) 
and, from the definition of W, symmetry under the interchange (p^,/^) «-> (p^., 
These represent further constraints on A, . . . , £ (eg. graviton indice symmetry requires 
(AjM+i = Ajk,s+2 & A, s +ij>+2 = Aj )S+2j j )S+ i) . However the Ward identities are 
sufficient to force A, . . . , E to zero even without such constraints and so these constraints 
are not explicitly incorporated into (|A.71|) for simplicity. We define 



^+i=^+2 = -E^ (A.72) 



and thus p$ +1 = p^ +2 give the momentum of the external graviton field. 

We also define Pj = Pi/iQ^Pi) for i — 1 . . . s + 2 and thus Q^Pt = 1 for all i. Then the 
homogeneous Ward identities corresponding to gauge and general coordinate invariance 
of condition (D) are given by 

Pt'W^,..^ = i = s + 2 (A.73) 

[Note that Weyl invariance implies a further homogeneous Ward identity that must be 
satisfied (which is not required in this section but will be used in the next subsection)] 

f +1 ^^..,, +lfc+J = (A.74) 



Now apply the Ward identities given by ( |A.73j ) to ( A.71 ). First examine the contribution 



which has (s-3) factors of in the contraction of P± with W 



[^42345?7/i2/« ? 7/-'4/i5 + ^2435^2/^4 

5^3/^4 n — - — - + 

[^2346^2/^3^4/^6 + ^2436^2^4^3/46 + ^2634^/i2/J6 7 7/'3M4] ( 5- " > + 

^1/^2/^3/^4/^6 

. . . + I [^s-l I s J s+l,s+2''7/i s _i/t a ''7/t i , + i/i s+ 2 "I" -^S— l,S+l,S,S+2 ? 7/i s _ip s+ i ? 7/ta/ts+2 

^.S-l,S + 2,S,S + l ? 7/t s _l/t s + 2 ? 7/i i ,/ia + l]Q' * -. | = 

^l^s-l^s^s+l^s+2 



However all terms above correspond to independent tensors and hence all the Ayki in 
the above expression must be zero. Hence all Aijki with none of {ijkl} equal to one must 
be zero. Contracting with the other P, shows that all A^i must be zero. 
Now consider the contribution with (s-2) factors of in the contraction of Pi with W . 
This gives the (with i,j not equal to one) in terms of linear combinations of the 
A^ki which from above are zero and hence such P^ must be zero too. Explictly the Pi 
contraction gives 



B 2 3V^2^Q- — " — > + B^^^Q, — * — , + . . . + Ps+i^+2^,+1^+2^ 



^1/^2/^3 ^1/^2/^4 ^+1^+2/^3 

+ I Terms involving Aijki = from above j = 

Considering all Pj contractions shows that all P^- are zero. By examining contributions 
with increasing numbers of it is straightforward to show that all (CV,-) Mfc , D Hllj and 
E^- are zero as well, which proves the desired result. Note that at 0(g 2 ), ie. the s=2term, 
the above method is invalidated because in order to show the Aijki are zero it is required 
to consider terms with (s-3) factors of Q^. 

We have now shown that the momentum space integrands of X HTL and X 9uess are the 
same at O(k). However, performing the momentum space integrals required to calculate 
X HTL and X guess is only well defined in Euclidean space for arbitrary fields. In Minkowski 
space the integrals are only defined for arbitrary fields via analytic continuation. However, 
for the fields belonging to 1Z, the support of the fields' Fourier transforms is sufficiently 
restricted for the Minkowski momentum space integrals to be well defined even without 
analytic continuation. Hence for fields in 1Z agreement of momentum space integrands 
gives equality of X HTL and X 9uess . [For fields not in 71, the same analytic continuation 
would have to be used for both X HTL and X guess to ensure equality]. 

A. 2 Relaxing condition (iii) above 

We also briefly remark that if condition (iii) above is relaxed then A^q,^ is not forced to 
zero by conditions (i) and (ii). The most general form of W^ 2fl3/14 taking into account 
the further constraints mentioned above is written down as W£p(P, R, K) as given by 
equation (Bl) of j7| noting the following notational changes 

([Pi,H ; [p 2)/ i 2 ] ; [p 3 ,h ; [p 4 ,H) <- ([R,t4 ; [-PM ; [K,a] ; [K,0\) 

Of course, this agrees with (|A.71|) for s=2 on imposing the additional constraints discussed 
above. Following the analysis of M we find that equation (B2) of M holds but equation 
(B3) of 0] is erroneous. Writing all the equations required to explicity show this is tedious 
and not very informative. However the end result , contrary to the erroneous conclusion 
of 0, is that the most general form of the 0(g 2 ) contribution to AS(C) a p consistent with 
the Ward identities is given by 

AS(C) a p = {gTf J dSdyd^re-^y-^y'-^A,{y)[C W^r)\A v {y') (A.75) 



with 



W$ = [{Q v P a - <)(Q^ - rft) + (Q"Pp - rjfii&Ra - rfc)] (A.76) 
-[Q{aV X p } ~ Q a QpK x ]KP[{Q»P p - V ;)(Q^R X - ttf) + (P, v) # (R, »)] 
~[Q a QpK 2 + Va p - Q {a Kp } ][(Q»P p - V $(Q"R>> - 

and C is either a constant or possibly a scalar dimensionless and rational function of 
the {(Q -p), (Q ■ r), (Q ■ k)} as multiplication by any other function will contradict (A) — > 
(D) of condition (i) above. 

Note that the coefficient of ( |2.32| ) can be written in the form of (|A.75Q for C = 
(Q ■ k) 2 /[(Q -p)(Q ■ r)]. This shows that the auxiliary field construction of section 2, which 
for E j£ satisfies (i) and (ii) but not (iii), is consistent with the above. The same applies 
for section 3. 



A.3 The Quark Sector 

In this subsection we outline a proof of the following: any expression, denoted 

X guess(quark)(Q. ^ ^ ^ ^ which 

i. Satisfies conditions (A) — > (D) 

ii. Reduces to X^ TL [A, t]] when e = 77 

must in fact equal to x HTL ( quarkS) (Q; tp, ip, A, e] at O(k) for all orders in g which is con- 
sistent with the assertions made in |J. We proceed in the same way as above and first of 
all define S 9 a J ss{quark) by 

s9 uess {q uark) = ( g a SX^^Q*,*^ (A>?7) 



Again, we suppose two different X guess (Q; A, e] exist satisfying (i)— ►(iii) above, leading to 
two different S^ ess whose difference we denote by AS^ ark \ The most general form of the 
0(g s ) contribution to ^s^ ark ^ can b e written as follows (with a — > fi s+ i and (3 — > /j s+2 ) 



A^ s+1Ats+2 = Jd A z 1 d 4 z 2 d 4 y 1 ...d 4 y s+2 d A r 1 d 4 r 2 d 4 p 1 ...d 4 p s (A.78) 



where jy s (<? uarfc ) i s a function of r 1; r 2 ,pi,p s+ 2 which are the momenta of, respectively, the 
anti-quark,quark and the s gluon fields. ji s+ \, jj, s+2 are the vierbien indices. The most 
general form of \y s =°(i uark ) consistent with conditions (A) — > (C) is given by 



s+2 



Hi 



+ E G^pIPifcTwT/* + EjkV^M (A.79) 

s+2>j>fc>l ^i/^fe 



where A' , A, Bj,Cj, Djk, Ejk are coefficients formed from the various contractions of 
the external momenta with the following possibilities: the external momenta, the gamma 
matrices or the <5 M .(Note that they contain no Q).The are proportional to some linear 
combination of the external momenta. The exact from of the above coefficients is heavily 
restricted by conditions (A) — > (C). 

Note that due to the relation 7^ 2 7 m = —7^7^ + 2?7 M1At2 t ne most general of form of 
yys=o{quark) can ^ e wr j^ en without terms proportional to 7 At .-7 Mfe with j < k and also with 
all $ dependence written with a single $ exclusively on the right hand side of QA.79Q . 

Now consider the homogeneous Ward identities. First we define 

vUi = Ps + 2 = -£-=irf - < - A (A.80) 

and hence p^ +1 = Pg +2 gives the momentum of the external vierbien field. Then with the 
definition Pf = (Q ■ for i — 1 . . . s + 2 the homogeneous Ward identities are given 

by 

P^W S =^;% =0 for j = 1 . . . s + 2 (A.81) 

H s+ lHs + 2-rxrs=0(guark) _ q jj^s(quark) _ jj^s(quark) /a oo 1 ) 

For the s = case, it is just an exercise in linear algebra to show that these identities 
force W s =^ quark) to zero. Thus we have the equality of X HTL( - quark ^ and X 9uess(quark ^ at 
0(g 2 n) (or rather equality of their momentum space integrands). 

The method of determining equality at higher orders in g is now identical to the 
method of the gluon sector i.e. consider the Ward identities given by QA.81 ) and the 



coefficients of the terms with increasing numbers of Q M factors. This is sufficient to force 

yys{quark) ^ ■ ^ crl UO n Case. 

Ati---A t s+2 ° 

Thus we have equality between the angular integrands of x HTL( - quark ^ and x 9uess ( quark ^ 
at 0{k) for all orders in g and thus X HTL{ - quark ^ = x gues < quark ^ at O(k). [See remarks at 
the end of section (A.l)]. 

B Appendix B 

The generalisation of the effective action to include the interaction of weakly coupling 
gravitons is independent (at O(k)) of the Weyl scaling properties of the auxiliary fields. 
Suppose, as in section 2, we have the auxiliary fields V(x, Q) and W(x, Q) transforming 
on Weyl scaling as — > e 2ua W^, and V — > e 2wT V^ with 
Then, the Weyl invariant extension becomes 



g 2 T 2 X E ^ v = J d 4 x(^g^ e A )2tr 



+ e~ vA V,y x F uX - V u [y a V a (e- uA W, + e~ uA {W u V ,y v - y a W a A„ )] 
+E{e-^ A V u [W x V,y x - y ■ WdA, + g, x W ■ Vy x + W,y ■ dA - y,W a A«}} 
dVL 



T E ee g 2 T 2 J^X E (Q;V,W,A,g] 



Notice that at O(k) the coefficient of E does not change as A ~ Vy ~ 0(k). 
Thus, we have 

X EjUjV (Q;V,W,A,g] = X EjU=0 ,v=o(Q; e~ vA V, e~ uA W, A, g] 

= X Efl>0 (Q;V,W,A,g])-K J dxkiuV ■ ^ + vW ■ J^) + 0(k 2 ) 

Note that the O(k) term on the right of the above vanishes at solutions of the auxiliary 
field Euler-Lagrange equations. Hence on eliminating the auxiliary fields by use of the 
Euler-Lagrange equations this term automatically drops out. Thus we have equality at 
(O(k)) between the two angular integrals on elimination of the auxiliary fields. This idea 
is easily applied to the other angular integrands of sections 3 and 4 as well. 



C Appendix C 

In this appendix we discuss further ambiguities that occur in attempting to write down 
T HTL [A, g] and Y HTL ^ quark > [A, e] to 0(k 2 ) using their known flatspace contributions and 
their symmetry properties given by (A) — > (D) above. Firstly consider A, defined to be a 
scalar which vanishes when g = i] with the Weyl scaling property that A —>■ A + 2a . 

This does not fix A uniquely at 0(k 2 ) and there are at least two degrees of freedom 
in A at this order. Consider A(oi, a 2 , x) defined by 



A(a 1? a 2 ,x) = -/ d6e Aiai ' a2 ' yix ' e)) de'e Aiai ' a2 ' yix ' e ' ) He- 2Mai ' a2 ' y{x ' e ' )) 



~RMx, &))y»{x, 6')f(x, 6') + i jfelM)) 2 + ai J 1 + a 2 J 2 

where J, = \V a y a {x, 9')V^(x, 9') - V a /(x, 6')V p y a (x, 9') 
■h = V^ a (x, 9')^y a (x, 9') - Vpy a (x, 9')\7 a f(x, 9') 

It is easily confirmed that that on Weyl scaling J\ and J 2 both merely scale by e 4<J 
and that f d9e A{ai > a2 ' y( - x ' 0)) is a Weyl invariant if A — > A + 2a on Weyl scaling [as a null 
geodesic, although invariant, is reparameterised on Weyl scaling with an affine parameter, 
9, transforming as d9 — > e~ 2a ^d9\. 

Thus, using the well known (see eg. Weyl scaling properties of R^, it can be 

shown that on Weyl scaling 

A(ai, a 2 , x) — > A(ai, a 2 , x) + 2a(x) (C.84) 

for arbitrary a\ and a 2 . Thus we see that there is at least a further 2-parameter ambiguity 
at 0(k 2 ) in trying to determine T HTL [A,g] and Y HTL ^ quark \A } e\. 

Note that A(0, 0, x) is the same as the expression used for A in || [as given by equa- 
tion^!) on using equation (55)] whereas A(l, 0, x) reduces to A = + J °° d9' V u y u \ y ( x ,e') 
ie. the A used throughout this paper. 

We also note that the local auxiliary field Lagrangians might well be dependent at 
0(k 2 ) on the Weyl scaling properties of the auxiliary fields suggesting the presence of 
further ambiguities at 0(k 2 ). 



D Appendix D 

Here we show that for any x 



[e-f I (y ■ V)" 1 ^]* = 



where we recall that 



with ^ vhc 
It is sufficient to show that 



(y-V)x = y v (d v + %gA v + uj' vhc o hc ) 



[e-~l (y-V)]x = 



(D.85) 



(D.86) 
(D.87) 



(D.88) 



(as (y ■ V) is invertible, we can write x = {v ' V) -1 ^ and then operating on ( |D.88j ) by 
(y ■ V) _1 e A gives (|D.85| ) for arbitrary 

Note that the following will be useful in going from the first to the second line of the 
calculation below. 

[ 7 a , J vhc o hc \ = +2^ a c7 c e aa rd,y a = -2/ Vr> aQ (D.89) 
where the right hand equation is due to the geodesic equation for y tl . Therefore 



[e~*j, (V ■ V)] X = [e-*y x e aXl \nd, + u' vbc a bc )]x 

= e"Vl/ A 7 c [r; A ea + e cX - e cX , ft + 2e aX u' u a c7 c ]x 

= e-t^Vh(e C A; M - 2e aA {+|e^e cp; J) 

+§(e cA A )jU -e CAt A, A +g fl xe?A, p )] 

Terms involving A drop out (note that g^xy^y = 0) and it is easy to check that the 
remaining term is identically zero. 

Hence [e"^, (jj ■ V)]x = as required. 
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